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ABSTRACT 

Thorne and Kidder give expressions which allow for analytical estimates of the 
"kick" , i. e. the recoil, produced from asymmetrical gravitational radiation during the in- 
teraction of black holes, or in fact any gravitating compact bodies. (The Thorne-Kidder 
formula uses momentum flux calculations based on the linearized General Relativity of 
gravitational radiation.) We specifically treat kicks arising in the binary interaction 
of equal mass black holes, when at least one of the black holes has significant spin, 
a. Such configurations can produce very large kicks in computational simulations. We 
consider both fly-by and quasicircular orbits. For fly-by orbits we find substantial kicks 
from those Thorne-Kidder terms which are linear in a. For the quasi-circular case, 
we consider in addition the nonlinear contribution {0{a'^)) to the ki cks, and provide 



a dyri amical explanation for such terms discovered and displayed by iBoyle &: Kesden 
(120071 ). However, in the cases of maximal ki ck velocities, the depend ence on spin is 



largely linear (reproduced in numerical results iHerrmann et al.l (l2007al )) 



Subject headings: black hole physics 
waves — relativity 



galaxies: nuclei — gravitation — gravitational 



1. Introduction 



Gravity waves (gravitational radiation) are a product of an extreme gravitational environment; 
the waves propagate through spacetime itself as fiuctuations in the gravitational field (the cur vature 
of spacetime). Experiments, supported by large collaborations such as LIGO (|Thorndll996l ). now 
seek to detect these waves of gravity by measuring these subtle fluctuations, further motivating 
theorists to understand gravity wave emission. Not only is it important to characterize these waves 
in the context of gravitational wave observations, but also for cosmolog i cal models in which black 
hole mergers play an increasingly important role (jdi Matteo et al.l 120081 : IWhitaker &: van Dokkum 
20o3). 



We consider two interacting black holes in a binary system that radiate away their energy 
and angular momentum with gravity waves. In addition to the orbital angular momentum of the 
binary, the individual black holes can also have their own spin angular momenta. If the black holes 
are gravitationally bound their radiation will eventually lead to inspiral and collapse. The final 
coalescence of t wo black holes r elease s a huge amount of gravitational radiation (up to 10 % of the 
total rest mass (jWashik et al.l (|2008l ))). If any sort of asymmetry is present in the binary, e.g. if 
the black holes have unequal masses, or if the black holes have unequal spin angular momenta or 
spin angular momenta unaligned with their orbital angular momentum, then the asymmetry will 
be reflected in the gravitational radiation emission of the coalescence, resulting in a "kick" of the 
final black hole. In this paper we consider the equal mass case of black hole merger where at least 
one black hole has substantial spin. 

Not only can we detect the gravitation radiation of the merger, the kick itself may propel the re- 
sultant black hole completely out of the center of a galaxy. Recent numerical simu lations of mergers 



of sp inning black holes result in kicks up to 4,00 km/s in quasicir cular inspiral (jCampanelli et al. 



20071 ) and 10,000km/s in hyperbolic encounters (IHealv et al.l 120081 ') . which could easily exceed the 
escape velocity of a galactic nucleus. IShields &: Bonnind (J2008l ) have presented limits on observing 
this phenomenon since kicked supermassive black holes may retain a portion of thei r in-falling mat 



ter, w hich may produce large flares of energy in a characteristic spectrum. Recently, iKomossa et al. 



(J2008l ) have discovered strong observational evidence of a recoiling supermassive black hole with 
optical emission lines. The black hole appears to have a kick of 2650 km/s! 



2. Multipole Formula 



Before the modern methods of numerical relativity were developed. iThornd (|l98Cl ) and iKidder 



(|l995l ) developed a multipole formula which describes the gravitational radiation kick from dynam- 
ical gravitating systems. The formula is based on derivatives of low-order multipoles of the masses 
and spins of the binary: 
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We have included only the terms that depend on the spin. Eq ([T|) is made up of nth time- 
derivatives '-"') of mass quadrupoles and octupoles, lij and lijk respectively, and of spin quadrupoles 
and octupoles, Hij and Hijk respectively. While the mass quadrupole and octupole are fairly 
familiar, the spin quadrupoles and octupole are less so. We present formulae below for these 
quantities, and a scheme to compute the spin multipoles. We use the physical Kerr parameter 
a = i /m, where j is the angular momentum and m is the mass of the black hole, and we point out 
that Hij and Hijk are linear in the spin a. Thus in Eq ([T]) the first and third terms are linear in 
the spin, and the second term is quadratic in the spin. 

Because Eq ([T]) involves time derivatives, it requires knowledge of the motion of the black holes. 



But there is no analytical form for two-body motion in General Relativity, so we will use Newtonian 
physics to describe the motion of the interacting black holes. (One effect of this choice is to confine 
the motion to the initial orbital plane, z = 0.) Thus our results will be accurate but uninteresting 
for Newtonian motions (slow velocities, large impact parameters or orbital separations). For the 
interesting case of relativistic interactions (velocities near c, small impact parameters or orbital 
separations) we will obtain qualitative results, which nonetheless lead to estimates of kick ratios 
for different configurations and allow an analytical under standing of the k ic k pro cess. For the 



quasicircular case our work is complementary to that of ISchnittman et al.l (|2008l ). which is an 
extensive study of non-equal mass mergers with either zero or equal but opposite spins; we study 
equal mass mergers with arbitrary ratios and directions of spins on the two black holes. 



2.1. Mass Quadrupole and Octupole 

The standard expressions for the mass quadrupoles and octupoles are: 

lij — / ^ Tn[XiXj —r OjjJ, 

^ijk = '^ rn[xiXjXk- -r'^i6ijXk + 5jkXi + 6kiXj)], (2) 

where the sum is over both black holes. We use Latin letters for spatial indices, x, y, z, and do not 
distinguish between covariant and contravariant indices. We treat one black hole at a time. The 
two-hole result is obtained by adding the contribution from each hole. Because in our examples the 
black holes have equal mass, the mass octupole vanishes on this addition, so that the third term in 
the formula Eq ([T|) can be ignored regardless of spin configuration. Since we need the third time 
derivative of the mass quadrupole, we evaluate Xk{t) and its first three time derivatives. 



2.2. Spin Quadrupole and Octupole 

The spin multipoles require the baffling concept of spin density, but since we are dealing with 
Kerr black holes, which have an intrinsic spin dipole moment, we use a trick to evaluate the spin 
quadrupole. We replace the spin dipole by a fictitious pair of spin cha rges, of value q = ifca sep arated 



by a distance m, centered at the actual location of the black hole (jHerrmann et al.ll2007al ). This 



reproduces the dipole angular momentum (am), and allows us to compute the quadrupole directly: 

Hij = ^q [xiXj - -r'^Sij] 
^ijk = X^ 9 [xiXjXk - -r'^{5ijXk + 5jkXi + 5kiXj)] . (3) 

Now the sum is over the two black holes and over the two spin charges for each hole, and the 
Xi appearing in these formula are offset in the direction of the spin. Strictly one should take a limit 



to small separations with the product am held constant, but this is automatic here because the 
result is proportional to am. We can separately consider the cases for spin components, ax, ay, a^, 
and add the spin quadrupole or octupole components after individual calculation. (To the lowest 
order, the spin is parallel-transported along the orbit.) 



3. Hyperbolic Fly-by 

Consider equal mass black holes approaching each other with impact parameter 26 in the x- 
y plane, and equal and opposite velocities (in the center of mass frame) of vq. We assume the 
motion starts along the x— axis, and the impact parameter is so large that the angle of deflection is 
small. Then, rather than explicitly computing the hyperbolic orbit, we can simplify the analysis by 
assuming that the motion is uniform in the x— direction. We will however compute the acceleration 
(and higher derivatives) using this uniform time dependence. Also, although the radiation reaction 
can change the orbital plane, we will not consider this feedback in the computation of the radiation. 
These points should become clear as we work through the analysis. In this section we consider only 
the lowest, linear- in-a term in Eq ([T]). 

Let us consider only the black hole initially moving along the line y = +b = const with velocity 
V = vq. It is at X = at time t = 0. The acceleration is given by the Newtonian result: 

(Px m,x m^VQt 

^~~4^~~4(62 + (^ot)2)|' 

where 2r = 2(6^ -|- {vqI)'^)^ is the separation between the black holes. Similarly, 

cPy m,h mh , , 
z. = = (5^ 

df^ 4r3 4(62 + (^ot)2)f 

(Our units have the Newtonian constant G = 1, and we henceforth take c = 1 also.) Integrating 
Eqs ([!])& dS]) once in time yields: 

Vx = vq^ -, (6) 

4T;o(62 + (^oi)2)2 



and 
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Similarly we can differentiate Eqs ([6]) & ([7|) to obtain higher time derivatives of the position. 



(7) 



In the Newtonian approximation the motion remains in the x — y plane (z = 0), and the motion 
of the equal mass black holes is symmetrical through the origin. Thus the contribution to the mass 
quadrupole is equal for the two masses. And, because of the symmetry, only I^x, Ixy, lyy, and I^z 
are in principle nonzero. In particular, Ixy = — ^ mxy and Izz = —Yl ■grnix'^ + y^)- 



Note that the deflection angle is of 0{m/bvQ). Consistent with our approximation, we keep 
only the lowest powers of {m/b) in computing multipoles. Our approach is encapsulated in the 
following rules: 

1. write the desired derivative of the multipole in terms of derivatives of x and y; 

2. replace undifferentiated x factors by vqI; 

3. replace undifferentiated y factors by 6; 

4. replace x (first time-derivative of x) factors by vq] 

5. in any term with a product of derivatives, first apply rules (3) and (4) above, then drop any 
term with more than one remaining differentiated factor. 



This approach is similar to those of lOohara &: Nakamural (|l989l ) and lBlanchet et al.1 (|l99Cll ). 



(3) 

We demonstrate the approach by evaluating Ixy using our prescription. We treat only one 
hole (the one with v^ ^ +vq): 
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As another example, we explicitly evaluate I- 
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Introducing the notation w 
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T-, the triply differentiated mass quadrupole for one black hole is: 
m'^vo w{2w'^ + 11) 
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(10) 



and other components are zero. These mass quadrupoles are for one black hole only, so in work below 
we include the contribution of the second equal mass black hole, which doubles these moments. 



We work out ^Hxx explicitly to demonstrate the method introduced in Section § 2.2. We 
assume spin on only one black hole. (To indicate the direction of the spin component generating 
the spin-quadrupole, Ux in this example case, we include a leading label ^ on the symbol ^Hij.) For 
the particle with velocity Vx ~ +^^0) we have: 

^Hxx = a.((x + f)'-^((x + f)2 + y2)) 
Aaxmx 



Similarly, 



'Hyy = a.(y2-i((x + |)2 + y2)) 

-ax{y'-l{{x-^f + y')), 

2axmx 
3 

X TT 

— J^zz, 
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= aziny, 



and others zero. 

After differentiating and combining components, we have: 

nxx - 3 l^flx^^g ay^^^} 



(11) 



'Hxy = ttxiny, (12) 

and others zero. (We used z = 0.) Also, 

yRxx = a.(x2-i(x2 + (y + ^)2)) 
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The linear term (the first term) in Eq ([1]) gives, for instance the force dPx/dt on the binary system: 

dPr 16 



dt 45 



(43)//£) + (43)-/f;)F(3)). (16) 



The force computed is apphed to the total mass, so we find the individual black hole velocity 
by time integrating Eq (fTUj) . using dt = b/vo, and dividing the result by 2m. The velocities are 
estimated in km/s for Ux ^ az ^ m, v ^ 1, and ^ ~ | (closest approach = 4m): 

t;^; ss 1962 fcm/s. (18) 



Similarly, 



^ = l^ff/(3) _J(3)^o-(3) _j(3)it{3)n 
J. /(rVv zz XX J xz xy yz /' 
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^ = l^f/(3) _ r{3)AiT(3) , r(3)^TT(3) _ it(3)n 
J. ^r v xx -^yy J^^xy '^ -^xy \ yy ^^xx )i 

_ ^ 2 m^vl b p ax{1w^+A) ^ ay{2w^ + hw) ^ ^ 

' ~ ^^ Ah b^ t;oloo (^2 + 1)4 + (^2 + 1)4 «^i 

= 2 X (-^(!p)V.o), (20) 

Vz ~ -4415 km/s. (21) 

Odd integrands integrate to zero in our straight-line integration approximation. The accumulated 
X— velocity is the residual CM motion after the encounter, but it is at most of order 10~^ of vq. 



so may be unmeasureable. However the Pz estimate of Eq ()2ip is substantial. Notice that these 
estimated kicks are for one black hole with spin. If both are spinning, the symmetries of the equal 
mass orbit dictate that Oj — > (oi — a2)i- Hence equal magnitude oppositely directed spin doubles 
this kick velocity. 

At this point we recall the limitations of these calculations, principally that the calculation 
of the dynamics is Newtonian. Our result is completely consistent and accurate in the Newtonian 
small-deflection limit, but the estimate Eq (j2ip is an extravagant extrapolation to vq = c. In the 
absence of a General relativistic 2-body simulation, we can make only qualitative adaptations to 
relativity. One point to notice is that -y is half the deflection angle in the high-speed Newtonian 
limit. For a test body moving near v = c past a central mass, in General Relativity the deflection 
at a given impact parameter and mass is twice the Newtonian result assuming v = c. This suggests 
that we might obtain the result estimated above from motion with twice the impact parameter. 



4. Quasi-Circular Inspiral 

To contrast our fly-by calculations above, we now calculate the kicks when equal mass black 
holes {mi = m2 = m) are in a circular orbit in the x-y plane. (In fact, the loss of energy means 
the orbit spirals inword, so is only quasi-circular, but we assume a circular orbit, with the orbital 
separation an adjustably shrinking quantity to mimic this energy loss.) We choose the first black 
hole initially (t=0) at position x = +d, where the second black hole would be at x = — d, with 2d 
as the "circular orbit separation" . The third derivatives of the mass quadrupole components for 
just the first black hole are thus: 

1^3) ^ -Amd^uj^sm.2Lot, 

/g) = -Amd^oj^ cos2ojt. (22) 

and other components are zero. The total differentiated mass quadrupole for the two equal-mass 
system is twice that given in Eqs ()22p . 

The spin multipoles are calculated with the method described in § 2.2 so that for one black 
hole, the two spin charges per component can be summed using the following coordinates: 

For Qt'- 
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xi = dcosujt, yi = dsuiijjt + {m/2), zi = 0, 

X2 = d cos tot, y2 = dsincut — {m/2), Z2 = 0. (24) 



For ttz 
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Again, as in the fly-by case above, we use the fact that the spin is to lowest order parallel 
transported along the orbit, which is (quasi-) circular here. The non-zero, third derivatives of the 
spin quadrupoles are thus: 



^hI^^J = -axdmuj^ siuLot, 
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and the non-zero, fourth derivatives of the spin octupoles, needed only for the second term of Eq 
(HD, are: 
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""H^J^J^ = Sazd^mcu'^ cos 2ujt, ^H^^'tl = -Sa^d'^muj^ cos 2Lot. (27) 



(3) 
Note that if both black holes have spin, ^H\' is computed by subtracting a similar formula for the 
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second spin: ^H^f^ ex (a^ — Ocj); and '^H^^L oc (ac^ + ac-j)- 



4.1. First Term 



We now calculate the first term from Eq ([T|), 



dP,. 



16 



dt p* 45 



««./f4?' 



(28) 



1* term only: 
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The first term of Eg dp is linear in spin, consistent with computational simulations as seen in 
Herrmann et alj (|2007al lbl) . The " 2 x " accounts for the two black holes of the system (doubling 
the mass quadrupole of a single black hole) . 

For arbitrary orientation of spin of magnitude a, the components are simply ax = a sin 9 cos (p, 
ay = a sin 6 simp and a^ = a cos 9. Above, the Oj are the spin components of just the one black hole 
that is spinning, but if both holes were spinning, we replace a with (oi — 02)- 

The circular orbit case presented here is based on Newtonian orbits, which specify frequency 
as a function of the Newtonian separation: 
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4.2. Second Term 

With the symmetries of lHerrmann et al.1 ( 2007a| ]bl). the second, nonlinear term in Eq ([T]) van- 
ishes identically. However if the spins are not equal in magnitude or not anti-aligned, this nonlinear 
term does not vanish, implying a quadratic contribution to kick velocity. 

We calculate the second, quadratic, term from Eq ([1]), 
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2nd ^gj.j^ only: 
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16 



= d m CO az{llayCosujt + 5ayCos3ujt — llaxSmivt + 5axSm3ujt). (32) 

The total force of the kick-components (first and second terms of Eq ([T]) together) can be 
compared in Fig ([1]) for the case where only one of two equal mass black holes is spinning at 
a = 0.6m with a separation of 6m. Notice that when the spin is perpendicular to the plane {9 = 0), 
the z-component of the kick vanishes. When the spin is in the plane {9 = ^), the linear term predicts 
only a z— component to the kick. In fact the in-plane kick components {x & y) are nonzero due to 
the non-linear contribution found in the second term of the multipole analysis formula (Eq ([I|)). 

With any choice of the orbital phase and in any spin configuration the largest kicks are those 
linear in spin. The contributions quadratic in spin become comparable to the linear terms only for 
a ~ 1. However we are interested in the general solutions rather than simply the largest kicks, and 
there exists a range of spin angles where the quadratic behavior of the second term in Eq ([1]) can be 
seen to dominate over the first, linear term in some components . This occurs when the spin angle 
9 satisfies |^ — fl < j^, that is, when the x-y spin components (orbital plane) are much larger than 
the z spin component. In this configuration the the x-y component of the kick is dominated by 
the quadratic contribution, while the z-component of the kick is more than an order or magnitude 
larger. 

Without an analytical 2-body solution, there is substantial ambiguity in converting this to a 
relativistic specification. An even more serious problem arises from the kick formulae predicting 
forces [dPi/dt) whose vertical component oscillates and whose in-plane components rotate with the 
orbit. If the orbit were strictly circular, the average of the kick would be zero, though the system 
would execute periodic motion due to the asymmetric radiation. In fact the orbit is only quasi- 
circular, and shrinks slowly due to gravitational radiation. Eventually the holes either spiral until 
they disappear behind a common black hole horizon, or enter a final plunge to the horizon. This 
sudden cutting off means the net kick can be modelled by considering the "last" (quasi-) circular 
fractional orbit. Such a concept is ambiguous at best because no analytic prescription describes 
the motion. However, numerical experiments do show kicks, and we can extract their dependence 
on orbital parameters. It is of great utility give simple analytical predictions, and thus we produce 
net, effective kicks depending on the phase of the orbit as it finally merges, and parametrized by 
the fraction of the "last" orbit that contributes the net kick. 
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Fig. 1. — Quasi-circular kick-components when the spins (a ~ 0.6m) are perpendicular to the plane 
{6 = 0) and in the plane {6 = ^). The x-component of the kick is the solid line, the y-component 
is the dotted line, and the z-component is the dashed line. Notice the quadratic effects appearing 
when the spin is in the plane. Quadratic effects will become comparable to the linear effects only 
when a ~ 1. 



12UU 






T 


1000 










- 








...■■E 




] 








,.E 


J 




- 


800 


- 


/ 


/'" 






- 


600 


: 1 








: 




/ 

T .■•■ 


/ 




' 




[ 


]■■ 




- 


400 








- 


700 




1 





6 (radians) 



Fig. 2 . — Comparing computational result for kick velocities (squares - the iHerrmann et al 



(l2007bll computat i onal B-Series) as a function of spin, to two analytic expressions: Dotted line 



Boyle fc KesdenI (120071 ) expansion; solid grey line - multipole analysis (this work). Because of our 



ambiguity in the "last orbit" radius, and in the phase of the fraction of the "last orbit" determining 
the kick, we fit two parameters: the maximum kick, and the phase. 
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Discussion 



It is of interest to compare the kicks from the quasicircular orbits, to the ones from the 
hyperbolic flyby. For instance, Eq(20) for hyperbohc flyby has a kick velocity 



4 m 
To compare this to the circular orbit result, consider the force, Eq.(29): 

dt ^6 ' 



(33) 



(34) 



where we used Eq(30). The kick velocity is approximated by multiplying -^ by a fraction (e/27r) 
of an orbital period, and dividing by 2m. Thus: 



la m d'^ 



I d^ 5 m^,a 



-^(-7) ( )Vorbit- 

22 d m 



(35) 



The process "multiplying -^ by a fraction of an orbital period" summarizes integrating the force 
for the relevant period at the final plunge, since circular orbits do not produce a net kick. 

Equations ()33p and (I35p are very similar, differing (aside from numerical factors) by vq in the 
hyperbolic case being replaced by evorbu, and the impact parameter by the orbital radius. The 
quantity e is poorly defined, but is likely less than unity. Also, while vq can in principle be very 
close to unity, the orbital velocity for a given orbital radius will be much less than the flyby velocity 
with an impact parameter equal to that orbital radius. This suggests that hyperboli c-orbit kicks 



can in principle be larger (much larger) than quasicircular kicks. Numerical studies (JHealv et al. 



(|2008l )) confirm these large kicks predicted by the multipole analysis for the fly- by case. 



Recently, iBovle Sz KesdenI (J2007l ) presented a spin expansion in order to understand final quan- 
tities of binary black hole mergers, such as mass, kick velocity, and spin vector. They consider two 
Kerr black holes in quasicircular orbits and Taylor expand some final quant ity in terms of the spins . 
Symmetry arguments remove excessive independent terms at each order. iBovle fc KesdenI (J2007l ) 
discover second and third order spin contributions that lie beyond the empirical fitting formulas 
which come from post-Newtonian, linear dependence fits from sii nulations. We coinpare o ur results 
to their expansion for a numeri cal black hole bin a ry sim ulation in lHerrmann et al.l (|2007bl ). referred 
to as the "B-series" (§IV C of iBoyle fc KesdenI (120071 )). In this particular case, equal mass black 
holes have oppositely directed equal-magnitude spins (a = 0.6m) lying in the x-z plane (6*1 =[0, 
vr] while 62 =[7r, 27r] for BHi and BH2, respectively). Fig ([2|) compares the Herrmann et al. nu- 
merically computed points with the Boyle & Kesden expansion and our multipole analysis for the 
resultant kick magnitude. Because of our ambiguity in the "last orbit" radius, and in the phase of 
the fraction of the "last orbit" determining the kick, we fit two parameters: the maximum kick, and 
the phase, and we plot the sum of Eq(29) and Eq(32). We find a tight agreement between these 
three methods, and the quadratic contribution in Boyle & Kesden expansion can be quantitatively 
understood with the quadrupole and octupole of the binary using the multipole formula. 
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